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)Q • Abstract 

C^, While presently available p-A data observe that the nuclear dependence of 

D , J/ip and ^p' production is the same within errors, we show that different 

nuclear dependence of J/^/J and V' production in the kinematic region un- 
covered by the present p-A data can be predicted based on two different 
nuclear absorption scenarios. It is found that the predicted production ratio 
a{tlj')/a{J/ip) at positive xp in A-p collisions is mainly determined by nuclear 
absorption, and hence allows direct experimental test of nuclear absorption 
scenarios. 
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Quarkonium production in experiments with a proton-beam incident on a nucleus tar- 
get, i.e. p-A collisions, has attracted much attention in recent years in connection with 
J/ip suppressions in high energy heavy-ion colhsions Actually it is a good place to 

check the hadronic absorption models which attribute the observed J/ip suppressions in 
nucleus- nucleus collisions to the absorption of J/ip through its interactions with hadronic 
environment instead of the formation of quark-gluon plasma (QGP). In particular, the exper- 
imental observations that the J/ ip and ip' production has the same nuclear dependence 
in p-A collisions put a strict constraint on various absorption models. However, it has been 
realized that all present experiments for the production of J/tp and ip' in p-A collisions 
studied only J/ip or ip' which is fast in the lab system or has positive Feynman xp. In 
fact, in the kinematic region not yet probed by present p-A experiments, i.e. the negative 
xp region, different nuclear dependence of J/ip and ip' production can be predicted based 
on different scenarios of absorption by target nucleons, which will be shown in this work. 
One of these scenarios, the "traditional" nuclear absorption model (we shall refer to it as 
^cenano-/ throughout this work), has been established in the literature [0-[T^ to account 
for the effects of absorption by nucleons in the following way: while the cc pair (pre- meson), 
produced as a colour singlet and of small spatial size, is expanding gradually to its bound 
state size, it can be dissociated by nucleons through an interaction cross section depending 
on the size of the cc pair. Very recently another scenario (Scenario-IT) stemming from 
the most recent theoretical progress made in hadronic production of quarkonium states, 
has been proposed to explain the observed charmonium suppression in hadron-nucleus and 
nucleus-nucleus collisions. This scenario considers charmonium production through the in- 
termediate next-to-leading Fock space state consisting of a colour octet cc pair and a soft 
gluon. The interaction between this composite state and nucleons can lead to charmonium 
suppression. In this work we shall show how these two different absorption scenarios could 
accommodate present p-A data and yet predict different nuclear dependence of J/ tp and ip' 
production in p-A collisions at negative xp- 

At negative Xi;' in p-A collisions the mechanisms of charmonium production such as 
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nuclear shadowing |T2|-|T4| and intrinsic charm p5|JT0|] are expected to be unimportant and 
will be neglected in our work. Initial state energy loss effects |T6| may lead to an increase 
of charmonium production in the negative xp region, which gives an opposite trend to the 



experimental data. Furthermore, the usual EMC effects ||T^ can be important at negative 



xp [|TB[. The initial state energy loss and EMC effects are not included in our study and their 
effects will be discussed later. Our previous work showed that Xi^-dependent comover 
contributions, as well as the nuclear absorptions, are essential for a comparison between 
theories and experimental data at small xp- We shall include only the effects of absorption 
by target nucleons and comovers in our study of J/ and ip' production in p-A collisions at 
negative xp. In fact, recent results on open charm production near xp=0 suggested that 
absorptions may be the dominant cause of the charmonium suppression in p-A collisions at 
negative xp. 

We shall show in this work that although the individual J/ip and ip' production at 
negative xp in p-A collisions predicted by the two different absorption scenarios may depend 
on the choice of parameters in comover contributions, the ratio of the two production cross 
sections (j{J /ip)/ a{ip') is insensitive to comover contributions, and therefore the prediction 
of the ratio may allow an experimental test of nuclear absorption scenarios. 

However, a study of charmonium production at negative xp (corresponding to the case 
of charmonium moving slow in the rest frame of target nucleus) is essentially impossible in 
an experiment of proton-beam incident on a nucleus target (p-A collisions) . An alternative 



is the experiment of a heavy ion beam incident on a hydrogen or deuterium target |[T8| , pT 
which we shall call A-p collisions. In A-p collisions the produced charmonia are moving fast 
in the lab system and thus the relevant measurements are more readily. In this work we have 
also calculated and compared the predictions of the mentioned two absorption scenarios for 
J/ip and ip' production in A-p collisions at CERN-SPS energy. 

In Scenario-I for p-A collisions the produced charmonium states {J /'ip or ip') or cc pairs 
(pre-mesons) can been dissociated by nucleons in the target nucleus into open-charm pairs 
(DD), which leads to suppressions of charmonium production. However, the formation of 



charmonium bound state requires finite proper time, and this formation time will suffer a 
Lorentz delay in the lab system. In the kinematic region reached by present experiments, 
the produced cc pair is moving fast in the lab system and may have been outside the target 
nucleus before it forms a physical charmonium bound state due to the Lorentz delay of the 
formation time in the lab system. Therefore, present experimental data show that the J/ if) 
and ip' productions have the same nuclear dependence in p-A collisions, although the radii 



of the J/ if) and i})' bound states differ by almost a factor of two in potential models 
According to this picture, those cc pairs which are slow in the lab system will have enough 
time to form the charmonium bound states {J/ip or ip' , etc.) inside the target nucleus, 
and thus the interaction of fully formed bound states J/ip and ip' with the surrounding 
target nucleons may lead to different nuclear dependence for J /ip and ip' production through 
the nuclear absorption depending on the final-state size, as we shall see below by explicit 
calculations. 

Including only the effects of absorption by nucleons and comovers, the a; i?- dependent 
cross section for quarkonium production in p-A collisions is 

^ ^ — oo I z 

xexp | - y ^ dracoVrenir, 6)| , (1) 

where a^^ is the bare nucleon-nucleon cross section, b the impact parameter and pA{b, z) the 
nuclear density. 

The second exponential factor in Eq. (|l]) represents the comover contributions, where 
(Tco is the (cc)-comover absorption cross section, Tq is the formation time of comovers, Tf 
is the effective proper time over which the comovers can interact with cc pair, Vre is the 
relative velocity of cc with the comovers, and n{T, b) is the density of comovers at the proper 
time T and impact parameter b. In our calculations we have introduced the xp dependence 



of comover contributions through n{T,b) as we did in Ref. |[T^, where n(T,b) is related 
to the Xi?- dependent comover distribution that is derived from the rapidity distribution of 
comovers. We fix Tf = ro/cs with rg ~1.2 fm and Cs — l/V^, fre —0.6 We shall vary 



the values of tq and (Tco to show the influence of comover contributions on the charmonium 
production at negative xp- 

The first exponential in Eq. (|1|) is the survival probability of the cc pair after its interac- 
tion with target nucleons. In the context of Scenario-I the nuclear absorption cross section 
is assumed to depend on the spatial separation of the cc pair, so that 

, , , <7RN{r/TRf if T < Tn, 

(^abs{z - Z) = < (2) 

aRN otherwise, 
where in general A; ^0 and we shall take fc=2 in our calculations. Here R stands for charmo- 



nium bound states J/ijj or ip'. The formation time of ^p' is taken as r^/ ^1.5 fm |2^, and the 
effective formation time of J/tp bound state is estimated to be ~1.2 fm, taking into ac- 
count the fact that about 70% of the observed J/ ip are directly produced and the remaining 



30% come from the decay of Xc |24]. We shall take the nuclear absorption cross sections for 
J /ip and Ip' as cr^Ar=4 mb from geometric considerations and a^i^ = cr^Nij^i /r^Y ~ Z.la^N 

©a- 

The proper time r in Eq. is related to the distance traveled by the cc pair through 
r = {z' — -2)/7f, V is the velocity of the cc pair in the rest frame of target nucleus, and 
7 = (1 — f ^)~^/^. We can relate the velocity of the cc pair in the rest frame of target nucleus 
to the Feynman Xi;' by 

P 



7^ = ^ = lc^{Pl + Pcm^iPD' + M^), (3) 

where P£ = XF{PL)max is the longitudinal momentum of the cc pair in the center-of-mass 
system and its maximum is 



{PDma. = ^[s - (M + 2m,Y][s - (M - 2m,Y]/2^s. (4) 

Here M is the mass of the cc system, irip the mass of a proton, and s the center-of-mass 
energy squared. The transformation from the rest frame of the target nucleus to the center- 
of-mass system is given by 7cm = (-E'b + "^p) / a/s, and (3cm = Pb/ (Eh + nip), where Ef, and Pb 
are the energy and momentum of the projectile beam in the lab system. One can see from 



Eq. that a. J/tp produced at negative xp corresponds to a slow J/ip in the rest frame of 
the target nucleus. 

We have calculated the xp dependence of a for J/ip and ip' production at negative xp 
in Fermilab E772/E789 p-A collisions in the parametrization, 

= ^^■"-). (6) 

dxp dxp 

The results are shown in Fig. 1 by the upper two curves, with the solid curve for J/%Ij 
production and the dashed for ip' production, parameters in comover contributions taken as 
crco=4 mb, ro=0.8 fm. Also shown are the Fermilab E789 data on the J/ip production at 
xp ^0 [Q. The curves are a little above the E789 data because for J/ip production the 



nuclear shadowing effect can contribute somehow at ~0 |]T9[. It can be seen that Scenario 



I actually gives the same nuclear dependence for J/ip and ip' production in p-A collisions 
aX xp ^0, as E772/E789 experiments have observed. However, different nuclear dependence 
for J/ip and ip' production appears at xp ^ —0.2, which is beyond the E772/E789 data. 
Note that the rise of aj/^ (solid curve in Fig. 1) at negative xp is due to the xj?- dependent 
comover contributions in Eq. (0). 

As we discussed above, in order to confront experiments we have calculated the produc- 
tion of J /ip and ip' in A-p collisions at CERN-SPS energy y^=17.4 GeV. As a matter of fact, 
the production at positive Xi;' in A-p collisions corresponds to the production at negative xp 
in p-A collisions. In Fig. 2 we show the results of J/ip and ip' production in A-p collisions 
at CERN-SPS energy for two sets of parameters cXco and Tq in comover contributions: the 
upper two curves with (Tco=2 mb and ro=1.0 fm (dot-dashed curve for aj/^ and dotted 
curve for a^'); the lower two curves with aco=^ mb and ro=0.8 fm (solid curve for aj/^ and 
dashed curve for a^i). One sees that comover contributions are essential for charmonium 
production at positive xp in A-p collisions. 

In order to reduce the uncertainties caused by the comover contributions, we have studied 
the production ratio R^pi = a{jp') / fa{J jip) in A-p colhsions at CENR-SPS energy, keeping 
in mind that comover contributions depends weakly on the type of produced charmonium. 



Here / ~ 0.14 is the experimental ip' and J/ip production ratio in hadron-hadron collisions 
p.|27[]. The results are illustrated in Fig. 3 by the upper two curves: the solid curve for 
parameters crco=4 mb, ro=0.8 fm, and the dashed curve for aco='2 mb, ro=1.0 fm. We see 
that the production ratio is insensitive to comover contributions. Moreover, in addition to 
nuclear absorption and comover contributions the EMC effects can be important at positive 



xp in A-p collisions ||T8[. However, since the EMC effects are initial-state effects and should 
be independent of the final produced type of charmonium, the inclusion of the EMC effects 
will not affect the value of the production ratio R^i The same argument exists for initial- 



state energy loss [|T^]. Therefore we suggest that the prediction of production ratio R^pi at 
positive Xi? in A-p collisions should be able to provide a direct test for nuclear absorption 
scenario discussed above. 

Another nuclear absorption scenario (Scenario-II) has been recently proposed by 
Kharzeev and Satz [0, and shown to be able to derive the phenomenological Gerschel- 



Hiifner fit |^ for hadron-nucleus data on J/ip production. In this scenario charmonium 
state is considered to be produced through the intermediate next-to-leading Fock space 
component | {cc)^g > formed by a colour octet cc pair and a gluon, whose lifetime is esti- 
mated to be Tg ~0.25 fm. The interaction of the composite state (ccg) with nucleons will 
lead to suppression of charmonium production. In particular, it is argued that the composite 
state {ccg) is of the same size for all charmonium states, and hence the (cc(7)-nucleon cross 
section is the same for all charmonia. As a result, the production of J/ip and ip' will have 
the same suppression, as long as charmonium suppression is dominated by the interaction 
of the intermediate (ccg) with nucleons. This is exactly what happens in the present p-A 
experiments. 

However, this scenario may also predict different nuclear dependence of J/ and tp' pro- 
duction in p-A collisions in the kinematic region where the composite state {ccg) has turned 
into the fully formed charmonium bound states J/ip and tp', and then nuclear absorption 
will depend on the specific type of charmonium state. This will be shown in the following 
calculations. 



Different from Scenario-I the nuclear absorption cross section in Scenario-II is now as- 
sumed to be 

(Tic-cg)N if Z'-Z< -fVTs, 

aabs{z - z) = < (6) 
o'RN otherwise, 

where jv is given by Eq. a^N are the relevant nuclear absorption cross sections as the 
same as those in Eq. (H), rg ^0.25 fm is the lifetime of the intermediate state (ccg) p9| , [Tl . 
The cross section of the interaction between the intermediate state (ccg) and nucleons cr(ccg)N 
is estimated to be 6 - 7 mb We shall take cr(^ccg)N=7 mb in our calculations. 



We have also in Scenario-II completed a calculation of the xp dependence of a in the 
parameterization Eq. for J/ip and ip' production in p-A collisions at Fermilab E772/E789 
energy. The results for chosen parameters, o"co=4 mb and ro=0.8 fm, are shown in Fig. 1 
by the lower two curves: dot-dashed curve for J/ip production and dotted curve for ip' 
production. One sees that at a;i;' ^0, i.e. in the kinematic region covered by E772 data on 
J/ip and ip' production, Scenario-II also gives nearly the same nuclear dependence of J/ip 
and ip' production. On the other hand, a rather different nuclear dependence of J/ ip and ip' 
production is predicted by Scenario-II at x^? ^ — 0.1. We can attribute this difference to the 
nuclear absorption dominated by the interaction of the fully formed J/ ip or ip' bound states 
with target nucleons. However, the predicted J/ip production at xp ^0 is more depleted 
than the E789 data, which might imply our over-estimate of the comover contributions in 
this region. The over-estimate, fortunately, will affect little the production ratio R^'/^, as 
it is shown below. 

In Fig. 3 we show the calculated production ratio R^p'/^ in A-p collisions at CERN-SPS 
energy based on Scenario-II by the lower two curves: dot-dashed curve with o"co=4 mb, 
ro=0.8 fm; and dotted curve with crco=2 mb, ro=1.0 fm. Again one find that the production 
ratio R^i given by Scenario-II is insensitive to comover contributions. Furthermore, the be- 
havior of R^' /.ij) given by Scenario-II appears quite different from that predicted by Scenario-I, 
which is due to the following essential differences between Scenario-I and Scenario-II: First, 
in order that the J/ip and ip' production be equally suppressed at positive xj? in p-A colli- 
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sions, Scenario-I permits very little nuclear absorption cross section at ^0 for both J/ip 
and ip' production, while in Scenario-II the interaction cross section between the intermedi- 
ate state {ccg) and nucleons is rather large but the same for J /ip and ip' production; Second, 
the proper formation time of physical J /ip and ip' , ^1.2 fm and t^i ~1.5 fm, is rather long 
compared with the proper lifetime of the intermediate state (ccg), rg ~0.25 fm, therefore the 
stronger suppression of ip' production relative to J/ip production in A-p collisions given by 
Scenario-II happens at lower value of xp than that given by Scenario-I (see Fig. 3). Based 
on the behavior of predicted R^'/^, an experimental test of nuclear absorption scenarios is 
possible. 

Although uncertainties of the parameters in nuclear absorption scenarios might reduce 
the power of prediction, some qualitative behaviors discussed here should be maintained by 
future experimental data. If not, e.g. the nuclear dependence of J/ip and ip' production 
remains the same at positive in A-p collisions, some nuclear absorption scenarios will be 
excluded as the explanation of quarkonium suppression observed in nuclear collisions. 

To summary, we have shown in this work that nuclear absorption scenarios are able to 
predict different nuclear dependence of J/ip and ip' production at negative xj? in p-A colli- 
sions or equivalently at positive Xi? in A-p collisions, although the present p-A data at xp ~0 
observed the same nuclear dependence of J/ ip and ip' production . Furthermore, we find 
that the production ratio a{ip')/a{J/ip) at positive xp in A-p collisions is mainly determined 
by nuclear absorption, but insensitive to comover contributions, initial-state energy loss and 
EMC effects, which can also contribute to charmonium production in this kinematic region. 
Therefore the predictions of production ratio at positive x^ in A-p collisions are suggested to 
be able to provide a rather direct and experimentally accessible test for nuclear absorption 
scenarios. 



9 



REFERENCES 

[1] T. Matsui and H. Satz, Phys. Lett. B178 (1986) 416. 

[2] M. C. Abreu et al. , Z. Phys. C38 (1988) 117. 

[3] J. Y. Grossiord et al. , Nucl. Phys. A498 (1989) 249c. 

[4] C. Baglin et al. , Phys. Lett. B220 (1989) 471; B255 (1991) 459; B262 (1991) 362; 
B270 (1991) 105. 

[5] D. M. Aide et al., Phys. Rev. Lett. 66 (1991) 133. 

[6] M. C. Abreu et al, Nucl. Phys. A566 (1994) 371c; C. Lourengo, Proc. Int. Workshop 
XXIII on Gross Properties of Nuclei and Nuclear Excitations, ed. H. Feldmeier and W. 
Norenberg (Gesellschaft fiir Schwerionenforschung mbH, Darmstadt, 1995) pl63. 

[7] S. J. Brodsky and A. H. Mueller, Phys. Lett. B206 (1988) 685. 

[8] J.-P. Blaizot and J.-Y. OUitrault, Phys. Lett. B217 (1989) 386; J.-P. Blaizot, Nucl. 
Phys. A498 (1989)273c. 

[9] S. Gavin and R. Vogt, Nucl. Phys. B345 (1990) 104. 

[10] R. Vogt, S. J. Brodsky and P. Hoyer, Nucl. Phys. B360 (1991) 67. 

[11] D. Kharzeev and H. Satz, preprint CERN-TH/95-214, BI-TP 95/30. 

[12] J. Qiu, Nucl. Phys. B291 (1987) 746. 

[13] J. Ashman et al, Phys. Lett. B202 (1988) 603. 

[14] S. J. Brodsky and H. J. Lu, Phys. Rev. Lett. 64 (1990) 1342. 

[15] S. J. Brodsky and P. Hoyer, Phys. Rev. Lett. 63 (1989) 1566. 

[16] S. Gavin and J. Milana, Phys. Rev. Lett. 68 (1992) 1834. 

[17] J.J. Aubert et al, Phys. Lett B123 (1983) 275; A. C. Benvenuti et al, Phys. Lett. 

10 



B189 (1987) 483; J. Ashman et al, Phys. Lett. B202 (1988) 603; P. Amaudruz et al, 
Z. Phys. C51 (1991) 387. 

[18] D. Kharzeev and H. Satz, preprint CERN-TH/95-73 BI-TP 95/16. 

[19] W. Q. Chao, C. S. Gao and Y. B. He, preprint CCAST-95-06. 

[20] M. J. Leitch et al, Phys. Rev. Lett. 72 (1994) 2542. 

[21] P. Hoyer, preprint NORDITA-95/78 P, |hep-ph/95114TI . 



[22] W. Kwong, J. L. Rosner and C. Quigg, Annu. Rev. NucL Part. Sci. 37 (1987) 325. 

[23] F. Karsch and R. Petronzio, Z. Phys. C37 (1988) 627. 

[24] S. Gavin, H. Satz, R. L. Thews and R. Vogt, Z. Phys. C61 (1994) 351. 

[25] B. Povh and J. Hiifner, Phys. Rev. Lett. 58 (1987) 1612. 

[26] M. J. Leitch et al, Nucl. Phys. A544 (1992) 197c. 

[27] R. Gavai et. al, in Hard Processes in Hadronic Interactions, Ed. H. Satz and X.-N. 
Wang; Preprint CERN-TH.7526/94, BLTP 63/94. 

[28] C. Gerschel and J. Hiifner, Z. Phys. C56 (1992) 171. 

[29] D. Kharzeev and H. Satz, Z. Phys. C60 (1993)389. 



11 



FIGURE CAPTIONS 



Fig. 1. Predictions for J /if) and if)' production in the Fermilab E772/E789 experiment 
with 800 GeV/c proton-beam incident on nuclear targets: the upper two curves given by 
Scenario-I (sohd curve for J /if) and dashed curve for ip') the lower two curves given by 
Scenario-II (dot-dashed curve for J /ip and dotted curve for ip'). Also shown are the Fermilab 
E789 data on the J/ip production p6| . 

Fig. 2. Predictions given by Scenario-I for J /ip and ip' production in the CERN-SPS ex- 
periment with 160 GeV/c P6-beam incident on a hydrogen target for two sets of parameters 
(Jco and To in comover contributions: the upper two curves with o"co=2 mb and ro=1.0 fm 
(dot-dashed curve for aj/^ and dotted curve for a^'); the lower two curves with aco=^ mb 
and ro=0.8 fm (dot-dashed curve for ajj^ and dotted curve for a^i). 

Fig. 3. Predictions of the production ratio for J /ip and ip' production in the CERN-SPS 
experiment with 160 GeV/c P6-beam incident on a hydrogen target given by Scenario-I 
(the upper two curves: the solid curve for parameters crco=4 mb, ro=0.8 fm, and the dashed 
curve for aco='i nib, ro=1.0 fm) and Scenario-II (the lower two curves: dot-dashed curve 
with aco='^ mb, ro=0.8 fm; and dotted curve with crco=2 mb, ro=1.0 fm). 
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